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Abstract. Present estimations of proton-proton total cross sections at very high en- 
ergies are obtained from cosmic rays (> 10 17 eV) by means of some approximations 
and the knowledge of the measured proton-air cross section at these energies. Besides, 
total cross sections are measured with present day high energy colliders up to nearly 2 
TeV in the center of mass (~ 10 15 eV in the laboratory). Here we use a phenomeno- 
logical model based on the Multiple-Diffraction approach to succesfully describe data 
at accelerator energies. Then we estimate with it proton-proton total cross sections 
at cosmic ray energies. On the basis of a forecasting regression analysis we determine 
confident error bands, analyzing the sensitivity of our predictions to the employed data 
for extrapolation. 



INTRODUCTION 



Recently a number of difficulties in uniting accelerator and cosmic ray values of 
hadronic cross-sections within the frame of the highest up-to-date data have been 
summarized [1]. Such united picture appears to be highly important for at least, 
the interpretation of results of new cosmic ray experiments, as the HiRes [2] and 
in designing proposals that are currently in progress, as the Auger Observatory [3] , 
as well as in designing detectors for future accelerators, as the CERN pp Large 
Hadron Collider (LHC). Although most of accelerator measurements of af% t and 
of^ at center of mass energy < 1.8 TeV are quite consistent among them, this 
is unfortunately not the case for cosmic ray experiments at y^s > 6 TeV where 
some disagreements exist among different experiments. This is also the case among 
different predictions from the extrapolation of accelerator data up to cosmic ray 
energies: whereas some works predict smaller values of af^ than those of cosmic ray 
experiments (e.g. [4,5]) other predictions agree at some specific energies with cos- 
mic ray results (e.g. [6,10]). Dispersion of cosmic ray results are mainly associated 
to the strong model-dependence of the relation between the basic hadron-hadron 
cross-section and the hadronic cross-section in air. The later determines the atten- 
uation lenght of hadrons in the atmosphere, which is usually measured in different 
ways, and depends strongly on the rate (k) of energy dissipation of the primary 
proton into the electromagnetic shower observed in the experiment: such a cascade 
is simulated by different Monte Carlo techniques implying additional discrepan- 
cies between different experiments. Furthermore, af^ in cosmic ray experiments 
is determined from cr l ^f: l air using a nucleon-nucleon scattering amplitude which is 
frequently in disagreement with most of accelerator data [1]. On the other hand, 
we dispose of parameterizations (purely theoretically, empirical or semi-empirical 
based) that fit pretty well the accelerator data. Most of them agree that at the 
energy of the LHC (14 TeV in the center of mass) or higher (extrapolations) the 
rise in energy of af% t will continue, though the predicted values differ from model 
to model. We claim that both the cosmic ray and parameterization approaches 
must complement each other in order to draw the best description of the hadronic 
cross-section behavior at ultra high energies. However, the present status is that 
due to the fact that interpolation of accelerator data is nicely obtained with most of 
parameterization models, it is expected that their extrapolation to higher energies 
be highly confident: as a matter of fact, parameterizations are usually based in a 
short number of fundamental parameters, in contrast with the difficulties found in 
deriving af% t from cosmic ray results [1]. If extrapolation from parameterization 
models is correct this would imply that o x ™ l air should be smaller, which would have 
important consequences for development of high energy cascades. 

With the aim of contributing to the understanding of this problem, in this paper 
we first briefly analyze in the first two sections the way estimations are done for 
proton-proton total cross sections from accelerators as well as from cosmic rays. 



We find serious discrepancies among both estimations. In the third section, on the 
basis of the Multiple Diffraction model applied to accelerator data, we predict pp 
total cross section values with smaller errors than with the standard techniques. 
We conclude with a discussion about the quality of present cosmic ray estimations. 

I HADRONIC aZ FROM ACCELERATORS 

Since the first results of the Intersecting Storage Rings(ISR) at CERN arrived 
in the 70s, it is a well established fact that af„ t rises with energy ( [11,12]). The 
CERN SppS Collider found this rising valid for of^ as well [14]. Later, the Tevatron 
confirmed that for af^ the rising still continues at 1.8 TeV, even if there is a 
disagreement among the diferent experiment values as for the exact value ( [15,16]). 
A thoroughful discussion on these problems may be found in [17,18]. It remains 
now to estimate the amount of rising of the total cross section at those energies. 
Let us resume the standard technique used by accelerator experimentalists [5]. 

Using a semi-empirical parameterization based on Regge theory and asymptotic 
theorems experimentalists have succesively described their data from the ISR to 
the SppS energies. It takes into account all the available data for a^t, P pp , &tot an d 
(? p , where pP p >P p , is the real part of the (pp,pp) forward elastic amplitude at t — 0. 
The fits are performed using the once-subtracted dispersion relations: 

C E r°° 
P±(E)a±(E) = ± + -/ dE'p' 

P Tip Jm 

where C s is the substraction constant. The expression for ailf v is: 

a% = AiE~ Nl ± A 2 E~ N * + C + C 2 [ln{s/s G )V (2) 

where - (+) stands for pp (pp) scattering. Cross sections are measured in mb 
and energy in GeV, E being the energy measured in the lab frame. The scale 
factor s have been arbitrarily chosen equal to 1 GeV 2 . The most interesting 
piece is the one controling the high-energy behaviour, given by a ln 2 (s) term, in 
order to be compatible, asymptotically, with the Froissart-Martin bound [19]. The 
parameterization assumes of^ and of^ to be the same asymptotically. This is 
justified from the very precise measurement of the p parameter at 546 GeV [28]. 
The eight free parameters are determined by a fit which minimizes the x 2 function 

The fit has proved its validity predicting from the ISR data (ranging from 23 to 63 
GeV in the center of mass), the of^ value later found at the SppS Collider, one 
order of magnitude higher in energy (546 GeV) [13,14]. With the same well-known 
technique and using the most recent results it is possible to get estimations for 
afo t at the LHC and higher energies. These estimations, together with our present 
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TABLE 1. a v t p ot data from high 
energy accelerators: fits values are 
from [5]. 
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experimental knowledge for both of^ and af^ are plotted in figure 1. We have also 
plotted the cosmic ray experimental data from AKENO (now AGASSA) [23] and 
the Fly's Eye experiment [24,25]. The curve is the result of the fit described in [5]. 
The increase in a^ t as the energy increases is clearly seen. Numerical predictions 
from this analysis are given in Table 1. It should be remarked that at the LHC 
energies and beyond the fitting results display relatively high error values, equal or 
bigger than 8 mb. We conclude that it is necessary to look for ways to reduce the 
errors and make the extrapolations more precise. 




II HADRONIC aZ FROM COSMIC RAYS 



Cosmic rays experiments give us af^ in an indirect way: we have to derive it 
from cosmic ray extensive air shower (EAS) data. But, as summarized in [1] and 
widely discussed in the literature, the determination of af^ is a rather complicated 
process with at least two well differentiated steps. In the first place, the primary 
interaction involved in EAS is proton-air; what it is determined through EAS is 
the p- inelastic cross section, of~™ r , through some measure of the attenuation of 
the rate of showers, A m , deep in the atmosphere: 

A - kX ■ — k 14 ' 5m P (A) 

I\ m — KA p - m r — K p _ air [4) 

The k factor parameterizes the rate at which the energy of the primary proton 
is dissipated into electromagnetic energy. A simulation with a full representation 
of the hadronic interactions in the cascade is needed to calculate it. This is done 
by means of Monte Carlo techniques [20-22]. Secondly, the connection between 
a ineT T an d &tot is model dependent. A theory for nuclei interactions must be used. 
Usually is Glauber's theory [7,8]. The whole procedure makes hard to get a general 
agreed value for af^. Depending on the particular assumptions made the values 
may oscillate by large amounts, from as low to 133 ± 10 mb [23] to nearly 165 ± 5 
mb [26] and even 1751$ [27] at ^ = 40 TeV. 

^From this analysis the conclusion is that cosmic-ray estimations of af„ t are not 
of much help to constrain extrapolations from accelerator energies [1]. Conversely 
we could ask if those extrapolations could not be used to constrain cosmic-rays 
estimations. 



Ill A MULTIPLE-DIFFRACTION APPROACH FOR aZ 

Let us tackle the mismatching between accelerator and cosmic ray estimations 
using the multiple-diffraction model [9]. The elastic hadronic scattering amplitude 
for the collision of two hadrons A and B is described as 

F(q, s)=i b db 1 - e* H(M) J (qb) (5) 
Jo 1 J 

where 5(6, s) is the eikonal, b the impact parameter, J Q the zero-order Bessel func- 
tion and q 2 = —t the four-momentum transfer squared. The eikonal can be ex- 
pressed at first order as 5(6, s) = (GaGb/), where Ga and G B are the hadronic 
form factors, / the averaged elementary amplitude among the constituent partons 
and the brackets denote the symmetrical two-dimensional Fourier transform. Given 
the elastic amplitude F(q,s), a\l t may be evaluated with the help of the optical 
theorem: 

<7^ = 4tt JmF(g = 0,s) (6) 



Multiple-diffraction models differ one from another by the particular choice of pa- 
rameterizations made for Ga and Gb and the elementary amplitude /. In the case 
of identical particles, as is our case, Ga = Gb = G. For our purposes we follow the 
parameterization developed in [10] which has the advantage of using a small set of 
free parameters, five in total: two of them (a 2 ,j3 2 ) associated with the form factor 
G 

G=U + £)-'<! + £)-' P> 

and three energy-dependent parameters (C, a, A) associated with the elementary 
complex amplitude / 

f(q,s) = Ref(q,s)+Umf(q,s) (8) 



1 - 31 



Imf(q, s) = C- ^ ; Ref(q, s) = X(s)Imf(q, s) 

1 _ „4 



(9) 



We get 



1 - e' n(b ' s) cos {\Q(b, s)}} b db J (q, b) (10) 



ImF(q — 0, s 

Jo 

with the opacity Q(b, s) given as: 

PCX) 

Q(b,s) = G 2 Imf(q,s) J a (q,b) qdq 

Jo 



(11) 



n(b, s) = C{E 1 K (ab) + E 2 K ([3b) + E 3 K ei (ab) + 

E 4 K er (ab) + b [EsK^ab) + E^K^b)} } (12) 

where k , k±, k ei , and k er are the modified Bessel functions, and E 1 to E 6 are func- 
tions of the free parameters. The proton-proton total cross-section is directly de- 
termined by the expression 

/CO f -\ 

b db[l -e^ (M cos[Afi(6,s)]} J (q,b) (13) 

This equation was numerically solved. The overall procedure is done in a three step 
process. 

• First, we determine the parameters of the model by fitting all the pp as well 
as pp accelerator data (differential elastic cross sections and p values), in the 
interval 13.8 < y/s < 1800 GeV. The obtained values are listed in Table 2. 



• Secondly, and most important, we estimate an error band for each of the 
energy-dependent parameters. To this end we introduce the so-called fore- 
casting technique, based on the multiple linear regression theory. It consists 
in determining a prediction equation por each free parameter. This allows 
to set a confidence band for each parameter and the confidence band for the 
predicted total cross section. The technique is explained in detail elsewhere 
[29]. 

• Finally, we proceed to extrapolate our results to high energies. Results are 
summarized in Table 3a and plotted in figure 2b, together with cosmic ray 
data. As a comparison, we list in Table 3b the extrapolated values obtained 
when only pp data, covering a much smaller t-range interval (13.8 < s/s < 62.5 
GeV), was used. This was the method in [10], but their extrapolated values 
were given without quoting any errors. 

It may be argued that af^t and af^ t are different at high energies: This is the "Odd- 
eron hypothesis", which as indicated in Section I, has been very much weakened 
[28]. Taking this into account, in our multiple-diffraction analysis it is assumed the 
same behaviour for af^ and of^ at high energy. 

Of course, if we limit our fitting calculations to the accelerator domain y/s < 62.5 
GeV (Table 3b), our results are the same as those obtained in [10]. In that case, the 
afo t values obtained when extrapolated to ultra high energies seem to confirm the 
highest quoted values of the cosmic ray experiments [26,27]. That would imply the 
extrapolation cherished by experimentalists is wrong. But the prediction af„ t = 91.6 
mb at the Fermilab Collider energy (1.8 TeV) is too high, first, and secondly, 
difficult to interpret, as no error is quoted in that work. In Table 1 we see that the 
measured of^ at 546 GeV is smaller than the predicted of^ by near 8 mb, and in the 
case of 1.8 TeV by more than 15 mb, which no available model is able to explain [18]. 
Also it can be noted that the extrapolation from figure 2a to ultra high energies 
may claim agreement with the analysis carried out in [26] and the experimental 
data of the Fly's Eye [27], and even with the Akeno collaboration [23], because its 
errors are so big that overlap with the errors reported in [27]. That is, such an 
extrapolation, Fig. 2a, produces an error band so large at cosmic ray energies that 
any cosmic ray results become compatible with results at accelerator energies, as 
it is claimed in [10]. However, when additional data at higher accelerator energies 
are included (Table 3b), both the predicted values and the error band obviously 
change. This can be clearly seen in figure 2b, where we have considered data at 
0.546 TeV and 1.8 TeV (see Table 1) in which case the predicted value of af^ t from 
our extrapolation at \fs = 40 TeV, o\l t = 131. 7t 4 ;g mb is incompatible with those 
in [26,27] by several standard deviations, though no so different to the Akeno results 
and the predicted value in [5]. 

Concerning the quoted error bands, the forecasting technique has reduced the 
errors, as is seen in figure 2(b), nearly by a factor of 3, as compared with the results 
quoted in Table 1. 



TABLE 2. Values of the parameters C, oT 2 , A at 
each energy. They are obtained from fitting the pp 
(pp) differential cross sections and p pp (p pp ) data in 
the interval 13 < ^ < 62.5 ( 546 < ^ < 1800 
GeV) 



V« (Tevj 


C (GeV ) 


a (GeV ) 
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13.8 


9.97 


2.092 


-0.094 


19.4 


10.05 


2.128 


0.024 


23.5 


10.25 


2.174 


0.025 


30.7 


10.37 


2.222 


0.053 


44.7 


10.89 


2.299 


0.079 


52.8 


11.15 


2.370 


0.099 


62.5 


11.50 


2.439 


0.121 


546 


??? 


??? 


??? 


630 


??? 


??? 


??? 


1800 


??? 


??? 


??? 











IV CONCLUSIONS: 

It has been shown in this work that highly confident predictions of high energy 
<j\l t values are strongly dependent on the energy range covered by experimental 
data and the available number of those data values. In particularly we show that 
if we limit our study of determining af^ at cosmic ray energies from extrapolation 
of accelerator data of \fs < 62.5 GeV, then results are compatible with most of 
cosmic ray experiments and other prediction models, because the predicted error 
band is so wide that covers their corresponding error bands (Fig. 2a). However, as 
x the included data in our calculations extends to higher energies, that is, when all 
experimental available data is taken into account, the estimated values for o\l t ob- 
tained from extrapolation and those obtained from cosmic ray experiments are only 
compatible, within the error bars, with the Akeno results (Fig. 2b). It should be 
noted that our predictions are compatibles with other prediction studies [5]. Taken 
all these convergences at face value, as indicating the most probable a^ t value, we 
conclude that if predictions from accelerator data are correct, hence, it should be 
of great help to normalize the corresponding values from cosmic ray experiments, 
as for instance by keeping the (k) parameter as a free one, as it is done for instance 
in [31]. The k value found will greatly help the tunning of the complicated Monte 
Carlo calculations used to evaluate the development of the showers induced by cos- 
mic rays in the upper atmosphere. In summary, extrapolations from accelerator 



TABLE 3. Predicted aZ from fitting accelerator data: 
(a) extrapolation including data at 546 GeV and 1.8 TeV 
(the two first values are interpolations), (b) extrapolation 
with data at \fs < 62.5 GeV; Experimental values are 
displayed in Table 1 
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FIGURE 2. Predictions (black squares) of oZ (a) width data at < 62.5 GeV; (b) including 
data at 546 GeV and 1.8 TeV. Open circles denote the interpolations. Notice the different vertical 
scales. 



data should be used to constraint cosmic ray estimations. 
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